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1 Introduction 



The knowledge of the hadron properties within the framework of Quantum Chromo- 
Dynamics (QCD) is one of the main issue in strong interaction physics. Several ap- 
proaches have been pursued in the last twenty years, in particular lattice gauge calcula- 
tions. Among the alternatives to these calculations, Light-Front Dynamics (LFD) is of 
particular interest []T|. It has proven successful in many phenomenological applications 
involving few-body systems in particle and nuclear physics. However, the application of 
LFD to field theoretical calculations is still in its infancy || . The main issue to be solved 
is the renormalization procedure ||. In perturbative calculations, the renormalization 
of the electron self-energy in QED, in standard LFD, is already non-trivial in the sense 
that it involves non-local counterterms |3J. This unpleasant feature is however a direct 
consequence of the choice of a preferential direction, the z axis, in the determination of 
the quantization plane. This can be well understood in the Covariant Formulation of 
Light-Front Dynamics (CLFD) ||, as shown in ref. 0. In this formulation, the state 
vector is defined on the light-front surface given by the equation uj-x = 0, where u is 
the four-vector with uo 2 = 0. The particular case where u = (1,0,0, —1) corresponds to 
standard LFD. In CLFD, the counterterm needed to renormalize the electron self energy 
in first order perturbation expansion is simply dependent on the orientation of the light 
front, defined by the four- vector 10. 

We shall investigate in this article how the question of non-perturbative renormal- 
ization can be formulated in CLFD. We shall first derive the general eigenvalue equation, 
whose solutions are the Fock state components. As a first example, we shall illustrate our 
strategy with a simple model involving two coupled scalar particles. Namely, the scalar 
"nucleon" N radiates the scalar "pions" ir. In this simple example, the Fock space is 
restricted to N, Nn and Nun states. Represented as a series of graphs in perturbation 
theory, it contains an infinite number of irreducible contributions to the self energy. They 
diverge and require renormalization. At large value of the coupling constant this system 
cannot be solved perturbatively. We show how to determine, in a self-consistent manner, 
the non-perturbative mass counterterm. This counterterm is then calculated numerically. 
Models involving two and three constituent bound states were also analyzed in 0. 

In a more general field-theoretical framework, in the Light-Front Tamm-Dancoff 
approximation involving spin 1/2 particles for instance, renormalization is not reduced 
to the introduction of a mass counter term. In this case, one should introduce sector- 
dependent counter terms, as shown in refs. ||. In the two-nucleon sector, for spin 
1/2 particles, additional box divergences appear ||. These divergences and the sector 
dependent counter terms are absent in the scalar model restricted to the dressed "one- 
scalar-nucleon" system. Therefore, their analysis is beyond the scope of the present 
paper. Inspite of that, the scalar model is still rather instructive, since the renormalization 
considered here is not reduced to the perturbative one. We will consider very large 
coupling constant, excluding the convergence of the perturbative series. 
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The plan of the article is as follows. In section ^]we establish the general equations of 
motion for the Fock components. In the truncated Fock space, the corresponding system 
of equations, which determine the Fock components and the mass renormalization, is 
detailed in sect. |3|. The renormalization of the wave function and of the coupling constant 
is calculated in sections £| and || respectively. Numerical results are presented and discussed 
in section [| We present our conclusions in section |7[ 

2 Eigenstate equation 

We start with the general eigenstate equation for the state vector B: 

P 2 <P(p)=M 2 <f>(p), (1) 

where 

h = H + p t nt ■ ( 2 ) 

We have decomposed here the momentum operator into two parts: the free one, P°, 
and the interacting one P*"*, given by: 

H = Z) / d t(k)di(k)k, d 3 k , 

i J 

f r + °° ~ dr 

P; nt = uj^ J H m \x)5{u-x)d A x = uJ ll J H mt (uT)—, (3) 

where we have denoted by H mt the Fourier transform of the interaction Hamiltonian: 

H int (p) = J H m \x)exp(-ip-x)d 4 x , (4) 

and d\ (dj) corresponds to the creation (destruction) operator for the various particles 
under consideration. The explicitly covariant formulation of LFD manifests itself in the 
fact that P^ n * in eq. ([3]) is proportional to uo^ and is determined by the integral over the 
light-front plane uj-x = 0. 

The equations for the Fock components can be obtained from (|l|), by substituting 
there the Fock decomposition for the state vector (j>{p) and calculating the matrix elements 
of P 2 in the Fock space. With the above expressions for P, eq. ([!]) obtains the form: 



(P ) 2 + (u>P°) / H m \ur)^ + I H^ioor^iu-P ) 



<j>(p) = M 2 4>{p) . (5) 
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In order to simplify this equation, we can use the fact that the operators (u-P°) and 
J H mt {ujT)dT commute. Indeed, from the commutation relation [P^, P v \ = we get: 

[u-P, P v ] = [u-P°, P° + F**\ = [u-P°, Pl nt ] = . 
Moreover, since uj 2 = 0, we can replace here and below (u>P°) by (uj-P). We thus obtain: 

(u-P°) [ H int (tur)(P(p)dr = [ {uj-P)H int {uT)(t){p)dr = (cu-p) [ H int {ujr)(j){p)dT , 



and the equation (|5|) is transformed to: 



dr 



2(u-p) j H m \ur)-<P{p) = - [(Py - M 2 \ 4>{p) . (6) 

The state vector <fi(p) is now decomposed in Fock components according to: 

d?h 



cj)(p) = (2n) 3 ' 2 Mh,P,uT)a\k 1 )\Q)6 iA \k 1 -p-uT)2(u-p)d 



T- 



(2^)3/2^/2^7 
+ (2ixf 2 J 2 (A: 1> A; 2 ,p > wr)at(fc 1 )6t(4)|O> 

x 8^ {h + k 2 -p- W r)2 M dr 

+ (2tt) 3 / 2 y fofa, h^^p^T^ihpfytffylO) 

x 5 {4) (k 1 + k2 + k 2 -p-ujT)2(uj-p)dr 
d 3 k\ d 3 k 2 d 3 k 3 

X (27T)3/2 ^ (27T)3/2 V ^ (27T)3/2^2^ + " ' ' U 



where = V 'hi + mf and is the mass of the particle i of momentum fcj. We introduce 
in (0) spinless particles of two types: the "nucleon" (with creation operator a') and the 
"pion" (with creation operator b>). The state vector (|7|) represents the "dressed nucleon", 
consisting of the "bare nucleon" and the admixture of one, two, . . . , many "pions". Our 
dressed nucleon is also scalar, therefore the state vector (|6]) corresponds to zero total 
angular momentum. This means that the Fock components are scalars and depend only 
on scalar products of all available four-vectors. 

The conservation law for the momenta in each Fock component has the form: 

k\ + k 2 H h k n = p + lut . 

Hence, the action of the operator (P ) 2 — M 2 in @ on the state vector <p{p) is reduced 
to the multiplication of each Fock component by the factor k^) 2 — M 2 = 2{uj-p)r. 
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If we introduce the notation: 



G{p) = 2(w-p)f0(p) , 

where f is the operator which, acting on a given component fa of 0, gives T0j. Q has 
therefore a Fock decomposition which differs from (|7j) by the replacement of the wave 
functions fa by the vertex parts Tj given by: 

r< = 2{uj-p)Tfa = (s- M 2 )fa , (8) 

where s = (XI ^i) 2 - We thus find the eigenvalue equation: 

i- 1 H m \ur)^g{ v ) = -g( P ) = -\(M 2 )g( P ) . (9) 

We introduce in (|9|) the factor A(M 2 ) depending on M 2 . The eigenvalue M 2 is found from 
the condition A(M 2 ) = 1. This equation is quite general and equivalent to the eigenstate 
equation ([!]). 



3 Equation for the Fock components 
3.1 System of coupled integral equation 

For the simplified model we consider in this study, we take the following interaction 
Hamiltonian: 

H = -g^ 2 (xMx) (10) 

where the scalar field ip with mass m corresponds to the scalar "nucleon" and the field p 
with mass fi corresponds to the scalar "pion". 

The system of equations for these vertex parts is shown graphically in fig. |l[ In 
order to write down this system of equations, it is enough to apply to the diagrams of fig. 
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Figure 1: Diagrammatical representation of the eigenvalue equation projected onto the 
Fock components of the state vector (restricted here to N=3). For the sake of clarity, we 
remove all kinematical variables. 



HI the rules of the graph technique detailed in ref. ||. We thus find: 

Y x {k x ,p^r x ) = 5m 2 / T 1 {k[,p,uT')5 {A \k[-p-uTy{u-k[)8{k'\-m 2 )d 4 k' l — 



t> -iO 



+ 9j^y 3 jv 2 {k' l X,P^')5^{k' 1 + k' 2 -p-ujT l ) 

x 6(LU-k[)5(k'l - m 2 )d%6(uj-k' 2 )8(k' 2 - ^ 2 )d%-^-^ , (11a) 



/j / 
r l (k[,p J ur')5^\k/ l -p-ury(uj-k[)5(k'l-m 2 )d 4 k[^- 

+ 5m 2 j T 2 (k[, k 2 , p, ut')5 {a) (k[ + k 2 -p- lot') 

x 6{u-k' x )8{k'\ - ™ 2 )d 4 k[-^-^ 



+ 



j T^k^k^p^r'^iK + h + k'z-p-LUT') 



x 9{uj-k' 1 )8{k'\ - m 2 )d%6(uj-k'z)5(k'l - ^ 2 )d%-^-^ , (lib) 



r 3 (fei,A;2, kz,p,UTz) = g J T 2 (k' 1 ,k 2 ,p,coT')8^(k' 1 + k 2 - p - oor') 



x eiu-k'Mk'l - m 2 )S%-^- 

T — IV 



Figure 2: Perturbative expansion, in terms of the pion-nucleon coupling constant g, of 
the nucleon self-energy. 



Figure 3: First order perturbative expansion of the nucleon self-energy. 



The origin of the mass counterterm Sm 2 in these equations is explained below. 
Since we truncate the Fock space to three particles, we omit in the last equation (|llc|) 
the coupling of to the four-body component Note that our approximation is not 
based on a perturbative expansion in terms of the coupling constant. It is based on 
a decomposition over intermediate states with an increasing number of particles. Once 
this number is fixed, we solve a non-perturbative problem in terms of g. If we iterate 
the system (|TTD, i.e., express T 3 and T 2 from the third and the second equations and 
substitute it in the first one, we get for Ti the sum of the diagrams shown on fig. 0. Of 
course, this series of graphs corresponds to all the irreducible contributions to the self 
energy with intermediate states up to Nnir. These contributions are iterated again, i.e., 
they appear repetitively on the nucleon line. In contrast to the case of the intermediate 
state Ntt, which generates only one self energy diagram shown on fig. |3], the number of 
irreducible contributions generated by the intermediate states up to Ntttt is infinite. The 
system of equations flTT|) corresponds to the sum of all of them. 

In the right-hand side of equations (|lla|) and ( lib ) we introduced the counterterm 
5m 2 corresponding to the interaction Hamiltonian H = —Sm 2 ip 2 (x). This term provides 
mass renormalization. The counterterm is however not introduced in the third equation 
containing the three-body intermediate state. At this point we have some freedom in 
the definition of the model Hamiltonian. In principle, the counterterm could be also in- 
troduced into the three-body intermediate state, since formally it does not increase the 
number of particles in the intermediate state. Our motivation in choosing the Hamilto- 



nian without the counterterm in the last n-th Fock sector is the following. Consider, for 
example, the mass operator —g 2 Y,(p 2 ) in second order perturbation theory. It is given 
by the diagram of fig. || with two particles (N and ir) in the intermediate state. The 
counterterm 5m 2 , which is shown by a cross on the nucleon line (one particle N in the 
intermediate state) just renormalizes this mass operator. So, one "pion" is deleted when 
the counterterm appears. Therefore, for the renormalization of the mass operator deter- 
mined by the sum of (an infinite number of) irreducible diagram with n particles in the 
intermediate state, one should consider the graphs, with the counterterm insertions, with 
n — 1 particles in the intermediate states only These diagrams are just generated by 
the Hamiltonian without the counterterm in the n-th Fock sector. 

We can easily transform eqs. ([□]) by performing the integrations which do not 
involve loops, keeping the loop integrals untouched. The result is the following: 

5 Tfl 2 

Ti(ki,p,UTi) = —. — 7— rx(ki,p,ujTi) 

+ 9j^y 3 Jv 2 (k[,k' 2 ,p,ujT')5^(k[ + k' 2 ~p-ujr') 

dr' 

x e(u-k[)S(k'l - m 2 )d%9(uj-k')5(k' 2 - fi 2 )d% , (12a) 

T — lO 



T 2 (k 1 ,k2,p,ujT 2 ) 
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+ 



2(uj-p)t 1 
5m 2 



Txip^p, W7i) 

r 2 (^i, k 2 ,p, ujt 2 ) 



2(u-p)t 2 x 1 



x 9(u-k[)5(k'l - m 2 )d 4 k[e(u-k' 3 )5(k' 3 ' - ^ 2 )d%- 
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dr' 



r'-zO 



(12b) 



r 3 (&i, k 2 , k 3 ,p,ujr 3 ) 



+ 



2(u-p)t 2 (1 - x 2 

9 

2(^)^(1-^3) 



-T 2 (k[,k 2 ,p, ujt 2 ) 
^2(k'(,k 3 ,p,UT 2 ) . 



(12c) 



In eqs. (12a) and ( 12b| ), 2{u-p)r\ = k\ — p 2 = m 2 — p 2 — > when p 2 — > m 2 . Since T\ is in 
the denominator, we keep p 2 7^ m? and take the limit p 2 — > m 2 in the final equation. In 
eq. (|12b|) , p\ = p — ojt\ with 2{u)-p)r 2 = Si 2 — m 2 and S\ 2 = (hi + k 2 ) 2 . In eq. (|12cj) we 
use the notations: 

2(uj-p)t 2 = s' 12 - m 2 , 2(lu-p)t 2 = s" 2 - m 2 , 

where 



'12 



(k[ + k 2 ) 2 , s'[ 2 



(k'i + h 
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r 2 . r\ . 
— C^y^ + — r ~*'' 

r 2 ^ r x _ ^ 

=^ = =^ ^ 



Figure 4: Diagrammatical representation of the eigenvalue equation in first order per- 
turbation theory. 

and k[, k" are determined by the conservation laws: 

k[ = k\ + k 3 + ojt' 2 - u;ri23, k" = ki + k 2 + ujt' 2 ' - ujt^ , (13) 

where ujt^ is the momentum of the spurion line entering the diagram. Everywhere we 
note Xi = uj-ki/uj-p . 

3.2 Reduction to two-particle Fock states 

We consider first the approximation in which the state vector (|7J) contains only the bare 
nucleon iV and the state Ntt. In this approximation, the system of equations for one- 
body and two-body Fock components is represented diagrammatically on fig. It is 
easily obtained by omitting T3 in ( pLlcj ), together with the counterterm in the equation 
which determines the last Fock sector T 2 . This truncation of the Fock space, retaining 
the minimal number of components, is equivalent to second order perturbation theory. 

Rewritten in terms of the variables q and n (see appendix A), this system of equations 
obtains the simple form: 

Ti = -4^F 1 + g f a(q',p 2 )r 2 (q',n) , (14a) 

T 2 (q,n) = — - Tx , (14b) 
m z — p z 

where cr(q,p 2 ) is the integrand which determines the self-energy XQo 2 ): 

/ • < 15 > 

1 1 

2 (S12 - p 2 - iO) e(q, m)e(q, ' 




S( P 2 ) = 
a(q,p 2 ) = 
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Figure 5: Same as Figure 1, but where the last Fock sector (n=3) has been expressed in 
terms of the two-body one. 

with 

S\2 = (h + k 2 f = [e(q, m) + e(q, /i)] 2 , e(q, m) = V 'm 2 + q 2 , (16) 

and similarly for e(q,fi). Since Ti does not depend on the relative momentum, it follows 
from eqs. that r 2 ((f, n) does not depend on the relative momentum too, i.e., r 2 (g, n) = 
const. Hence, we get: 

f-i + 4^Vi + ^( P 2 )r 2 = o, 

\ m z — p J 

/ 2 rx - r 2 = o. (17) 

From ([17|) we find the eigenvalue equation: 

p 2 = m 2 - 5ml - g 2 Y>(p 2 ) . (18) 

The counter term Sm^ is determined from the on-shell condition p 2 = m 2 , where m is the 
physical mass of the nucleon. This gives: 

5ml = -s 2 £(m 2 ) , (19) 

as expected in second order calculation of mass renormalization. 

3.3 Solution for the three particle system 

Consider now the system of equations which incorporates the three-body Fock component 
It was shown graphically on fig. 1. We can easily express T3 through r 2 . The system 
of equations is thus transformed as shown in fig. 5. Its analytical representation has the 
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form: 



r 2 (g,n) 



Ti+^y a(g',p 2 )r 2 (g',n) 

g ^ , 5m 2 + ^ 2 S(si) 

— r 2 (g,n) 



m 2 — p 2 



(2^)3 ' 



(20a) 



m 2 — p 2 



(Sl2 - P 2 )^fc! 



3 / n(g,g / ,n,p 2 )or(g / ,p 2 )r 2 (g / ,ra) 



(2tt)* 



(20b) 



Each term in eqs. ( PU| ) is represented by a graph on fig. 5. In these equations, s is given 
by ([16D, and: 

si = (Aii - wr) 2 = m 2 - x fcl (s i2 - p 2 ) , 



uj-p 



[e(q, m) — n-q\ 



In this equation, H(q, q',n, p 2 ) is the propagator of the second intermediate state: 



n(g,g',n,p 2 



dr" 



1 6 [u-{h - k' 2 )\ 5 Uh -k' 2 + lut 2 - ut") 2 - m 2 ] — - . 

J T — iO 

9 [uj{k[ - k 2 )} 



m 2 — (k[ — k 2 — lot') 2 



(21) 



Its expression in terms of the variables q, q', n can be easily calculated with the kinematics 
detailed in the appendix. 

The system of equation (EOT) can now be solved by two independent methods. 



i) The first one consists in the elimination of Ti in the two equations in (20). The 
result is an equation for T 2 which, at p 2 = m 2 , obtains the form: 



T 2 {q,n) 



g 



5m 2 



d 3 q' 
(2*0 



,P ) T 2(q ,n) 



+ 



5f 2 E(si) +5m 



[Sl2 ~PV x k 
2 



2 f 

—T 2 (q,n)+g 2 J U(q, q' ,n,p 2 )a(q' ,p 2 )T 2 (q' ,n) 



d 3 q' 
(2^)3 



= X(5m )T 2 (q,n) 



[22) 



Like in eq. (|Sp, we introduce in ( |2"2"D the factor A (5m 2 ). The mass renormalization 5m 2 
will be found from the non-perturbative condition X(5m 2 ) = 1. ii) The second method 
consists in the direct elimination of T 2 in the two equations (|20D . We can indeed rewrite 
eq. (|20b|) in the following form: 



n) — 



bra 2 + g 2 Y>(s 

(Si 2 ~p 2 )x k 



—T 2 (q,n) - g 2 J U(q,q\n,p 2 )a(q\p 2 )T 2 (q\n) 



d 3 q' 
(2^)3 
9 



m 2 — p 2 



ri . (23) 
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After discretization of the momenta, this equation can be written schematically: 

A l3 Y\ = ~^-^Y 1 [if , (24) 
m 2 — p 2 

where [1] is the vector whose every component is 1, and is a two-dimensional matrix 
obtained after discretization of the momenta (index % for q,n-q and index j for q',n.q'). 
The vertex function T 2 can now be expressed in terms of Ti after a simple matrix inversion: 

n = -jr-jiW . (25) 

m z — p z J 

After insertion in eg. ( |20a|) , we end up with an equation involving Ti only. Since T\ is a 
non-zero constant, it can be removed from the equation, leading to an equation for Sm 2 



generalizing eq. (|19|). In the perturbative limit, the matrix A reduces to the unit matrix, 



and we therefore recover eq. (|i~9|) exactly. 

4 Renormalization of the wave function 
4.1 Non-perturbative case 

We schematically rewrite the state vector \p), given by eq. (0), as: 

| J 9) = 1 |iV)+0 2 |iV7r) + 03|iVvr7r} . 
It is normalized as follows ||: 

(p'\p) = 2p 5^(p-p') . (26) 
The Fock components are thus normalized in order to provide the condition (26) . Substi- 



tuting the state vector (|7|) in the left-hand side of eq. fl2H|), we get, before normalization 
of the Fock components: 

(p'\p) = Z2po5^(p-p') , (27) 
Z = N 1 + N 2 + A<3, 

N x = <j>\, (28a) 



where 
with: 



k2 ,^ Je(q,m) +e(q,fi)] 3 



2(2tt) 3 J e(q,m)e(q,ii) 
N 3 = J ^I(qi,q2,q3,n)5 (3) (q 1 + q 2 + q 3 ) 

[e(q 1 ,m)+e(q 2 ,n)+s(q 3 ,fi)\ 3 3 3 



e(qj.,m)e(q 2 ,ii)e(q 3) fjL) 
12 



d i q 1 d i q 2 d i q 3 . (28c) 



From T 2 calculated with eq. (j2"2"D, we can find 02 according to: 

/_ »\ Tofo, n) 
S12 - m 2 

The calculation of iV 2 is then straightforward. To calculate N\ and N 3 we should know 0i 
and 03 . The calculation of 0i and 3 is explained in appendix [B]. Note that 0i does not 
depend on any relative momentum. The normalized state vector satisfying the condition 
(^6|) obtains the form: 

|p> = <P? n \N) + ^Mr) + fr\Nw) , 

where 

C£s = 4>x,2, 3 /Vz . (30) 

We can introduce the creation operator of the new, composed field, creating directly the 
state \p). It is written schematically as: 

A\p) = <^ e V + (jf^aW + 0™ n aW , (31) 

so that \p) = A(j9)|0 >. With the renormalized wave functions ren , the vacuum expec- 
tation value of the commutator: (0| [A(p), A^p)] |0) = 2p S^(p — p') is the same as the 
one-body operators a, a* (except for the normalization factor 2p ). So the state \p) which 
is a dressed state in terms of the bare operators ,b\ can be interpreted, in this sense, 
as an elementary particle in terms of the operator A*. 



4.2 Perturbative case 

The perturbative case is simply obtained from the preceding one by omitting 3 . Tak- 
ing into account that from eq. ([14]) we have Ti/(m 2 — p 2 ) = 0i, we find: T 2 = g<f>\. 
Substituting it into fl2"9"|), we find 02 and then, by eq. ([28 b|), we obtain N 2 : 



^ = 0^/2, where J 2 = / ^ ^ (32) 

Wir 6 J (su - m z Y e(q,m)s{q, /i) 

We thus find the ratio: 

^ = 167rm 2 a/ 2 « 0.38a (33) 

2 

with a = 16 ^. m 2 ■ The numerical value of J 2 is given for m = 0.94, fi = 0.14. The 
integral I2 is logarithmically divergent for \i — > 0. The two-body contribution is rapidly 
decreasing, when the mass of the intermediate particle increases. 
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5 Renormalization of the coupling constant 



As mentioned above, the mass renormalization counterterm 5m 2 , for the state vector in- 
corporating N and Ntt states only is given by eq. (|i"9|) and coincides with the perturbative 
result. For the state vector incorporating the states N, Ntt and Ntttt, it is determined by 



eq. (|22|) . This renormalization constant is infinite when the cutoff tends to infinity. 

The coupling constant is also renormalized, though its renormalization is finite for 
the particular scalar system we are interested in the present study. We show below how 
this renormalization can be carried out. This renormalization is a "by-product" of the 
renormalization of the wave function fulfilled in the previous section. Let us consider first 
the case of the state vector containing N and Ntt states only. 

5.1 Truncation to iV and Ntt states 

In this approximation, the normalization factor N 2 is given by eq. ([32|). We thus get: 

Z = N 1 + N 2 = (fiZt with Z x = 1 + g 2 I 2 , (34) 

where I 2 is given by eq. (^). After renormalization the wave function (p 2 turns therefore 
into: 



9<f>l j_ren 9 9r 



(35) 



s — m A \/Zi(s — m 2 ) s — m z 

where we introduced the renormalized coupling constant 

9ren = 9 1 (36) 

This value of g ren can be also represented as the residue of the two-body wave function at 
s = m 2 , i.e. the value of r 2 (g, z) at the nonphysical value of q corresponding to s = m 2 . 

One can alternatively define g ren from the irN scattering amplitude determined by 
N exchange in the s-channel: 



m? — p 2 m? — p 2 m 2 — p 2 m? — p 2 — 5m 2 — g 2 Tj(p 2 



where p 2 = (k\ + k 2 ) 2 — 2{uj-p)r. Near the pole p 2 = m 2 we get: 

77 9 9ren 



(m 2 — p 2 )(l + g 2 dY>(p 2 ) / dp 2 \ 2=m 2) m 2 — p 2 
where we introduced g ren by eq. (|36|) , but with Z\ given by: 

dnv 2 ) 

^1 — ^ + 9 — yi - 

ap z 



(37) 



(3? 



(39) 
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Taking eqs. flTSp for £ one can easily check that the equations ( pip and ( |3"9"D determine 
the same Zx. 

Note that the renormalized coupling constant g ren is finite and it is always smaller 
than the bare one g. When the bare constant g tends to infinity, the constant g ren remains 
however finite, but it reaches its maximal value. Expressed in terms of a = g 2 /(167im 2 ) 
it has the form: 

cCT = l/(16vrm 2 / 2 ). 

It does not depend on g. According to (|33|) , its numerical value for m = 0.94 and fi = 0.14 
is: cCf = 1/0.38 = 2.63. 

5.2 Truncation to N, Nir and Ntttt states 

In this case the renormalized wave functions are given by eq. (|30|) with Z determined by 
eqs. (^7],^). Note that the renormalization constant Z can still be represented in the 
form ( |3"9"| ) with S(p 2 ) determined by all the irreducible contributions (see appendix [C]). 

Restricting ourselves to the system considered above and having found the renor- 
malized state vector \p), we can get all the physical information (we can for example 
calculate the electromagnetic form factors, if the particles are charged). So we do not 
need in practice to define and calculate the renormalized coupling constant g ren . How- 
ever, it is useful to calculate it for further generalization to the case of particles with spin, 
when the charge renormalization constant will become infinite. 

As we already mentioned, the standard definition of the coupling constant is the 
residue of the wave function 2 at s = m 2 , i.e., the value of T 2 (q, z) at the non-physical 
value of q = in corresponding to s = m 2 . For a non-relativistic bound state calculation 
for instance, it is given by k = where is the binding energy of the bound 

state. In coordinate space it is the coefficient of the asymptotical behavior of the wave 
function ip{r — > oo) oc exp(— kt). Since we calculate T 2 (q,z) numerically in the physical 
region for q, it is not easy to find numerically with enough accuracy its interpolation into 
the non-physical region. However, we do not need to choose the renormalization point 
s = m 2 , but we can choose any other renormalization point. We can define for instance 
the renormalized coupling constant as the value of T 2 (q, z) at q = 0: 

~g 2 ren = Tr(q = 0,z). (40) 

Note that T 2 (q = 0, z) does not depend on z. In the case of iV + Nn intermediate states 
T 2 does not depend on q and the renormalized coupling constants g ren and g ren coincide 
with each other. 

6 Numerical results 

Both methods to solve eqs. (|TT|) are used to cross-check our results. For regularization 
purposes, we introduce a cutoff L, i.e. integrate in (|22|) over modulae of all the relative 
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three-momenta q until q < L. Note that this cut-off procedure preserves rotational 
invariance. For the masses, we choose the nucleon and pion masses: m = 0.94 GeV 
and \x = 0.14 GeV. The integration over the azimuthal angle is done analytically. The 
equation is reduced to a matrix form by discretizing the integral. The convergence of the 
integrals is already obtained for 30 points in the variable q and 15 points in the variable 
z. The points in q where not taken equally spaced, but with a spacing proportional to 
h 2 , where h is the equal spacing in the variable Jq. In the first method, the eigenvalue 
X(Sm 2 ) of the matrix is found numerically and 6m 2 is fixed to get A = 1 to less than 0.5%. 
For the second method, 5m 2 is calculated by a standard iteration procedure, starting from 
the perturbative result. Since both methods give identical results to less than 1% we only 
quote the results obtained with the second method. 

The results for the dimensionless coupling constant a = g 2 /(167rm 2 ) = 3 and for 
different values of the cutoff parameter L are shown in table |l|. We denote by 8m 2 /5ml 
the ratio of the value 5m 2 found from eq. (P2"[), to the perturbative value 5ml, given by 



eq. (|T9T). The values Aq j2j 3 denote the contributions of the corresponding Fock sectors 
to the normalization of the state vector. In order to characterize the wave function 
quantitatively, we calculate also the average value of the (absolute) relative momentum 
< q >, normalized to the two-body Fock component: 

N 2 2[2ny J e{q,m)e{q,n) 



L 


1 


5 


10 


50 


100 


200 


5ml 


-1.234 


-3.67 


-4.82 


-7.54 


-8.71 


-9.88 


5m 2 /8ml 


1.095 


1.052 


1.040 


1.025 


1.022 


1.020 


Nt 


0.377 


0.329 


0.327 


0.326 


0.326 


0.326 


N 2 


0.441 


0.455 


0.456 


0.458 


0.458 


0.458 


N % /N x 


1.183 


1.380 


1.395 


1.407 


1.405 


1.405 


N 3 


0.182 


0.216 


0.217 


0.216 


0.216 


0.216 


< q > 


0.366 


0.538 


0.565 


0.588 


0.591 


0.593 



Table 1: Numerical results for a = 3, as a function of the cut-off L. 

One clearly sees that though 5m 2 increases logarithmically, the contributions of the 
Fock components N±, N2, N3, as well as the average momentum < q > in the two-body 
Fock component become stable after L = 5. This means that we indeed found numerically 
the renormalized solution for the wave function. 

One can see also that the nonperturbative value of 5m 2 is very close to the perturba- 
tive one Stuq. We have checked that the function ^((7, z) is almost constant, as expected 
for the perturbative solution. The solution for Sm 2 remains very close to the perturbative 
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result also for higher values of the coupling constant. The results for L = 200 and for 
different a between 1 and 1000 are shown in table |2[ One can see that 5m 2 is very close to 
the perturbative value even for extremely large coupling constant a = 1000. As expected, 
the three-body sector dominates, when the coupling constant increases. However, the 
ratio N 2 /Ni is still close to its perturbative value (A^/iV 1 ) , given by eq. (|33l). 



a 


1 


3 


10 


30 


100 


1000 


5m 2 /5ml 


1.011 


1.019 


1.029 


1.035 


1.040 


1.043 


N, 


0.661 


0.325 


0.081 


0.016 


2. E-3 


3. E-5 


N 2 


0.292 


0.457 


0.366 


0.186 


0.068 


8. E-3 


N 3 


0.047 


0.218 


0.553 


0.798 


0.930 


0.992 


N 2 /N 1 


0.441 


1.41 


4.51 


11.91 


32.9 


230 


(Wo 


0.38 


1.14 


3.8 


11.4 


38 


380 



Table 2: Same as Table 1, but with L fixed to 200, and a is varied. 

The reason that the solution is close to the perturbative one lies in the super- 
renormalizability of the scalar theory. The third term in eq. (|22|) converges and does not 
require any cutoff, whereas the first two terms are divergent and dominate. Therefore, a 
very small departure of 5m 2 from its perturbative value 5ml can accommodate the finite 
higher order correction. If the third term can be neglected, the equation is approximately 
satisfied with the perturbative value of 5m 2 . In order to show that it is indeed so, we 
introduce the cutoff in the two divergent terms, but do not introduce it in the third, 
convergent one. The results for a = 3 are shown in table |3|. For small enough value of L the 
two first terms are suppressed, and the solution drastically differs from the perturbative 
on. When L increases, the solution becomes closer and closer to the perturbative one. 



L 


0.1 


1 


10 


100 


5m 2 /5ml 


36.3 


1.55 


1.07 


1.03 



Table 3: Test calculation in which the cutoff L is introduced only in the two first terms 
in eq. |2^ (see text), for a = 3. 

We emphasize that the small deviation of the solution from the perturbative one, 
even for very large values of the coupling constant, is just a property of the non- 
perturbative equation fl22"D , which includes the contributions (with one-, two- and three- 
body intermediate states) to all orders of g. This result cannot be justified in any pertur- 
bative expansion in terms of g. 
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7 Conclusion 



In a first attempt to address the question of non-perturbative renormalization in CLFD, 
we have investigated in this study a simple, but nevertheless meaningful, model based on 
two scalar particles. This model is reminiscent of the structure of the physical nucleon in 
the low energy regime, in terms of bare nucleons coupled to pions. 

Using the nice features of CLFD, we have first derived the general eigenstate equa- 
tion that should be used in order to calculate any physical state vector. We emphasize 
here that this equation is quite general and is not restricted to the case of scalar particles 
nor to the restricted Fock space we consider in this study. It should therefore also be used 
when solving more complex systems as QED or QCD. 

The results we obtained for the simple scalar model, with a restricted Fock space 
expansion up to three particles, are quite encouraging. We obtained, numerically, a renor- 
malized solution using a simple mass counter term. Surprisingly enough, this counter term 
is not very different from the perturbative, logarithmically divergent, mass counterterm, 
even for very large values of the coupling constant. We traced back this feature to the na- 
ture of our scalar model we start from, which is super-renormalizable. This result however 
does not imply that higher Fock states are negligible. We find that the first non-trivial 
Fock component gets larger and larger as the coupling constant increases. 

The direct generalization of this study is the investigation of non-perturbative renor- 
malization in scalar QED, following the study of perturbative renormalization in QED 
already done in ref. || in CLFD. This will be the subject of a forthcoming publication. 
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A Kinematics 

We show in this appendix, how to express the propagator (pip in terms of the variables 
q, q", n. 

The two-body wave function depends on the following four-vectors: 

(j) — <j){kx, k 2 ,p,ojr), k 1 + k 2 =p + uJT. (42) 
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We introduce the variables: 



n = L- X (p)i2/\L- X (p)i2\ = VV^L- x {V)uo/u-p , (44) 

where 

V = p + ujt , (45) 

and L^iV) is the Lorentz boost. The wave function under the integral (see the last 
diagram in fig. 3) depends on k[, k' 2 ,p,UT', and, correspondingly, on the variable: 

q' = L-\V')k\ , (46) 

where 

V' =p + ujt'. (47) 

In order to obtain II (<f, q', n,p 2 ), we should express the four-momenta k[, k 2 , ujt' in (k[ — 
k 2 — ujt') 2 in terms of q, q', n. Since under the Lorentz transformations and the rotations 
of the four-vectors the variables q, q', n are rotated only, the expression for the scalar 
II(g, q', n, p 2 ) does not depend on the system of reference, and, hence, can be found in 
the most convenient one. We find it in the system where 

p' = k' 1 + k' 2 =p + ujT' = . 

— * _ — * 

In this system k[ — q ' and k' 1Q = e(q',m). From the conservation law the vector V 
is expressed as: V = k\ + k 2 = uj{t — r'). We have: 

r = t' = (48) 

2(u- P y 2(uj-p) ' 1 1 



where yfs = e(q,m) + e(q,fi), vs 7 = e(q',m) + e(q',fi). Since uj = uj n and uj-p = 
uj-(k[ + k' 2 ) = uj (e(q', m) + e(q', pi)) = ujq\/~s' \ we find 

^ ^ s — s ' ^ s + s' , , 

V = n^^, Vo = • 49 



One can check that V 2 = s. According to eqs. ( |i3|) and (§§), the variables q and q' are 
defined by different Lorentz boosts. From ( fi"9"l) one can see that this difference is the boost 
in the direction of n. This boost does not change the unit vector n. 
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Now let us find k 2 . It is obtained by reverting (031): 



k 2 = L(V)(-q) = -q + - J=[e(g,/i) - -=J^ ] 

Vv 2 Vv^ + v 



(50) 



Eq. floTf ) is obtained from ( H3J ) by replacing in the r.h.s k\ by — q and by changing the 
sign at V . 

Substituting here eqs. (EST) for V, we find: 



-q + n 



L s — s 



2Vss 



e(q,fi) -n-q 



(51) 



and similarly for k 20 : 



e{q,fj,)Vo q-V .s + s' ^^s-s' 



2Vss 



2Vss 



From (|48| ) we get: 



cut = n ^=-, lu t 



, s' — p 2 



(52) 



2v^ ' 2^ ■ 

Substituting the above expressions for the four-momenta into (k[ — k 2 — ujt') 2 , we find 
that the four-vector squared in the denominator of the propagator ( pT|) is expressed in 
terms of the variables q, q', n as follows: 



(k[ -k 2 - UJT 1 ) 



>\2 



,->, \ / x S + s' _ _S — s' s' - p 21 
e{q , m) - e{q, fi) — p= + n-q- 



2 V / ^ 7 2v^s 7 2V 7 ? 



,s — s 



Zi \J s s 



— n- 



, s' — p 2 
2V 7 ? 



-I 2 



.(53) 



B Calculation of and 

The Fock components <pi and 03 are calculated as follows. The wave function <p 2 in 
(|28b ) is related to the vertex functions T 2 by eq. (p9| ). Having found numerically the 



vertex function r 2 (<f, n) = r 2 (g, z) from the equation (|2"2"D, we then find by fl2T3| ) the Fock 
component (f) 2 . 

With eqs. ( |12a| ), reduced to the first equation in (|14|) , we express Ti through Y 2 in 
the limit p 2 — ■> m 2 , and then find the component 0i: 



5m 2 



| ^{q\p 2 )Y 2 (q\n)^ . (54) 
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With eqs. ( |12cj ), we express T 3 through T 2 and then find <p 3 : 

, - ^ 9<t>2{q',z') g<f) 2 (q",z") 

(si23 - m 2 )(l - x 2 ) (S123 -m 2 )(l -x 3 ) 

Here <f) 2 (q',z') is <f>2(k' 1 ,k2,p,ojr£) represented through the relative momentum, and 

<h(q",z") is MK,fo,P,ur2). 

The wave function 2 in the r.h.s. of eq. ( [55] ) depends on the variables defined in 
the center of mass of the two-body subsystem. We shall now express these variables in 
terms of the three-body relative momenta q\, q 2 , (&■ 

The variable q' is related to s' l2 by s' 12 = [e(q',m) + e(q',jj,)} 2 . We thus get: 

g' 2 = [ m 4 + (/i 2 _ s , i2) 2 _ 2m 2 (fl 2 + s , 2) ] /(4s / 2) 

where, by taking the first of eqs. ([13]) squared: 

, S13 - m 2 

Su — s i23 : > 

1 - x 2 

S123 = (e(q 1 ,m)+e(q 2 ,n)+e( y q 3 ,n)) 2 , 
S13 = (fa + k 3 ) 2 = (e(q u m) + e(q 3 , /J,)) 2 - q\ , 

_ e(gi,m) -ra-gj _ e(g 2 ,3, v) ~ n-g 2 , 3 . , 

^1 — -j= j ^2,3 — -j= • I^OJ 

VS123 ySl23 

We use the fact that q\ + q 2 + q 3 = 0. 

The value of z' is found by comparing expressions (|5Tj| ) for xi with: 

efg', m) — z'g' 



with: 



£(q',m)+£(q',n) ' 
This gives: 

z' = [e(q', m) - x x [e(q f , m) + e(q', mu))] / q' , 

where X\ is given by fl5"E|) in terms of qi,q 2 ,q 3 . 

The values of g", 2" are found from the equations: 

q" 2 = [m 4 + (/i 2 - s'[ 2 ) 2 - 2m 2 (/i 2 + s'{ 2 )] /{As'[ 2 ) , 

z" = [E(q",m)-x 1 (e(q",m)+e(q",n))}/q" , 

where: 

// _ S12 - m 2 

s n — s i23 ; > 

1 - x 3 

with S12 = (A;i + /c 2 ) 2 = [e(qx,m) + e(q 2 , /i)] 2 - <? 2 . 
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C Proof of eq. ( ggj ) in the general case 



We show below that the renormalization constant Z can still be represented in the form 
(P9|) with S(p 2 ) determined by all the irreducible contributions. In this case, the number 
of irreducible diagrams containing two- and three-body intermediate states is infinite. 
Some of them are shown in fig. |[ Consider, for example, a contribution containing n 
successive intermediate states. The corresponding amplitude contains the factor 1/Tj for 
any of these states. Like in eq. (fL5l ), any factor 1/n turns into l/(sj —p 2 ), where Sj is the 
invariant energy of a given intermediate state and p is the external incoming/outgoing 
momentum. So, indicating only these factors, we represent this contribution to X(p 2 ) as: 



p2 



where the dots include all the integrations with the corresponding measures. Derivative 
over p 2 gives the factor l/(si — p 2 ) 2 - 

Now consider the graphs which are the same at the right to a given z-th intermediate 
state and differ from each other by the contributions at the left to this state. The infinite 
sum of them determines the amplitude of the virtual transition from the initial state N to 
the states Ntt or Ntttt. It is the vertex T 2 (if the state i is the two-body state), or T 3 (if 
the state i is the three-body state). We can then take the sum over all the contributions 
to the right of this given state i and again obtain r 2 or IV So, the result has the form 

*w_r n ... i .1 ... (.*....(*... v.. • a.. 



and, after extraction of the common factor Ni, we recover eq. fl39] 
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